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■ bstract. In this paper, we examine a generalized magma equation for rational values of 

two parameters, m and n. Firstly, the similarity reductions are found using the Lie group 
method of infinitesimal transformations. The Painleve ODE test is then applied to the 
travelling wave reduction, and the pairs of m and n which pass the test arc identified. These 
particular pairs arc further subjected to the ODE test on their other symmetry reductions. 
Only two cases remain which pass the ODE test for all such symmetry reductions and these 
are completely integrable. The case when m = 0, n = — 1 is related to the Hirota-Satsuma 
equation and for m = i, n = — i, it is a real, generalized, pumped Maxwell-Bloch equation. 
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1 Introduction 



O 
O 

Scott and Stevenson (22] examine the process of melt migration where buoyant magma rises 
through the Earth's mantle. The system is treated as a two-phase flow, with the porous rock 
matrix able to deform by creep as the melt rises. The model can be written in terms of the 
porosity or liquid volume fraction eft, i.e. volume of melt in unit volume of the fluid-rock mixture. 
The constitutive relations k = kQ<j) n and r\ = rjQ(p~ m are assumed, where k is the matrix per- 
r> | meability, rj is an effective viscosity for the matrix and fen and r/n are constants. After rescaling 4> 

with its uniform background level, so that the background now corresponds to (f> = 1, Scott and 



Stevenson [Jfy derive the partial differential equation 

A = 4> t + (<T + m^- 1 ^ - r~ m <Pt x ) x = 0, (1.1) 
and suggest that physically appropriate ranges for m and n are 

< m < 1, 2<n<5. 

The same equation has also been derived independently for m = by McKenzie ^5] . 

Equation is well known to have solitary wave solutions and various authors have exami- 
ned these. Nakayama and Mason |2J looked for rarefactive solitary waves for which > 1 
everywhere, by using a travelling wave reduction. This leads to 

° (£)'='<*>■ 
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for some function / where z = x — ct and ip(z) = <fi(x,t). They then applied the boundary 
conditions 

/(1) = 0, (1.2a) 

|(D=0, (1.2b) 

/W = 0, *>1, (1.2c) 

and proved that such waves exist if the index n > 1 and do not exist if < n < 1, with their 
analysis restricted to n > 0. They go on to derive some solutions for particular values of pairs 
(m, n) and to look at the large amplitude approximation. In contrast, in a later paper |23| . 
a perturbation solution for small amplitude rarefactive waves is examined and at leading or- 
der has the sech-squared form of a single-soliton solution of the Korteweg-de Vries equation. 
Nakayama and Mason |2U1 122j also look at compressive solitary waves, for which < < 1. 
For boundary conditions, they take the first two to be the same as those for rarefactive waves 
()1.2a|) . (|1.2b|) . but in place of the last, they instead use 

0(1) =0. (1.2d) 

Zabusky |34| and Jeffrey and Kakutani used this condition for examining algebraic com- 
pressive solitary wave solutions of the modified Korteweg-de Vries equation. However, they also 
found exponential solitary waves and it therefore appears that this choice of boundary condition 
is limiting. Indeed, in this paper we exhibit some new exponential solitary wave solutions. 

In Nakayama and Mason [20], with the boundary conditions (jl.2a|) . (|1.2b|) . (|1.2d|) . the wave 
speed c must be equal to the index n. They derive three compressive solitary wave solutions 
which tend algebraically to the background ip = 1 on either side and reach tp = at their 
maximum amplitude. The wave with n = m = 2 has a monotone behaviour whereas the one 
with n = m = | is oscillatory. The solution for n = 3 and m = is not differentiable at 
the point of maximum amplitude, where there is a cusp. In |22| they examine the necessary 
conditions for these compressive solitary waves to exist with the same boundary conditions as 
before. They find that it is necessary to have n = m > 1 and plot graphs of the curves for 
integer values 2, 3, 4 and 5 for m and n. It is found that the shape does not greatly change, but 
that the width increases slowly as n = m increases. The half-integer values |, |, |, and | are 
also examined and all prove to have similar oscillatory structure. 

Takahashi and Satsuma |29j make a change of variables and find a periodic wave solution for 
the choice n = 3, m = in terms of elliptic integrals. They also investigate weak solutions and 
demonstrate the possibility of both a hump solution and two travelling wave solutions between 
different levels with a sharp wavefront. They note that if the if the 4>t x term in (|1.1[) is replaced 
by 4> xx , it is possible to apply a transformation to yield the Korteweg-de Vries (KdV) equation. 
However, the transformation is not necessarily single- valued and the significance of this result 
is unclear. Takahashi, Sachs and Satsuma |2H] show that the travelling wave solutions with 
compact support suggested in [221 are n °t physically feasible, due to singularities in the stress. 
They use them, however, as suitable initial conditions for various numerical simulations and 
demonstrate the break-up into a series of solitary waves. 

Marchant and Smyth ^S] also examine periodic wave solutions for n = 3, m = and 
develop a modulation theory for slowly varying travelling waves. It is found that full or partial 
undular bores are possible and there is good comparison between approximate wave envelopes 
and numerical solutions. 

Experiments by Scott, Stevenson and Whitehead [2Zj) Olson and Christensen |24j and White- 
head [SJ show solitary waves on conduits of buoyant fluid in a more viscous fluid. The relevant 
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equation is identical to the magma migration equations for m = 1 and n = 2 and the experi- 
mental results are summarized in Whitehead and Helfrich [33]. They show Korteweg-de Vries 
'soliton-like' collisions (cf. JO]), in which two solitary waves interact and then emerge unaltered 
apart from a phase-shift. The experiments also imply that an arbitrary initial condition quickly 
breaks down into an ordered sequence of solitary waves. These results raise the question about 
whether the partial differential equation (jl.lj) is completely integrable (i.e. solvable by the in- 
verse scattering transform) for any of the parameters m and n. This is further reinforced by 
Whitehead and Helfrich [32] who demonstrate that (jl.lj) reduces to the KdV equation for small 
disturbances to the uniform state. 

Barcilon and Richter [3J have examined rarefactive waves for m = and n = 3 both ana- 
lytically and numerically. They concluded that this case did not have soliton solutions as there 
was an interaction left behind by the collision of two solitary waves and they could only find 
two conservation laws. A completely integrable partial differential equation such as KdV has 
an infinite number of such laws. 

Harris also examines conservation laws for (jl.lj) . without limiting the values of m and n 
to the physically relevant range, and proved that there were at least three laws for m = n + 1, 
n 7^ 0; at least two laws for m = l, n/0 and precisely two laws for all other combinations of m 
and n. Thus the special cases m = n + 1 and m = 1 (with n / for both situations) are the 
only ones which could possibly lead to completely integrable equations. 

In this paper, we apply Painleve analysis to look for any such completely integrable partial 
differential equations, omitting only the case m = n = for which the equation is linear. We 
start by using the Lie group method of infinitesimal transformations to find all the similarity 
reductions of the partial differential equation to ordinary differential equations. We can then 
apply the Painleve ODE test [2j to these and they must all have the Painleve property (no 
movable singular points apart from poles) for the partial differential equation to be completely 
integrable. 

Firstly, we examine the travelling wave reduction. Applying the Painleve ODE test, we 
find that there are five rational pairs (m, n) which pass for any value of the wave velocity c. In 
addition, there are also four cases which require c = 1, but these can be discounted as candidates 
for complete integrability. 

The five remaining cases are then subjected to the Painleve ODE test on their other similarity 
reductions. Three of these fail the test at this stage, but the other two pass the procedure. These 
last two cases are then checked with the Painleve PDE test [30j and are found to be completely 
integrable. 

One case, m = 0, n = — 1 is related to the Hirota-Satsuma equation (or Shallow Water Wave 

1 Equation) examined by Clarkson and Mansfield The second pair, m = |, n = — ^, yields 
a partial differential equation which can be transformed to a special case of the real, generalized, 
pumped Maxwell-Bloch equation. This has been previously investigated by Clarkson, Mansfield 
and Milne [7j . We note that both pairs satisfy the criterion of m = n + 1 given by Harris 

for the possibility of more than three conservation laws. 

2 Symmetry analysis 

The classical method for finding symmetry reductions of partial differential equations is the Lie 
group method of infinitesimal transformations. As this method is entirely algorithmic, though 
often both tedious and virtually unmanageable manually, symbolic manipulation programs have 
been developed to aid the calculations. An excellent survey of the different packages available 
and a description of their strengths and applications is given by Hereman j!2M13j . In this paper 
we use the MACSYMA package symmgrp .max [5] to calculate the determining equations. 
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Table 1. Infinitesimals for the magma equation Ql.lfl . 



m. 


n 






T 


77 


n ^ m. 




aix + ao 


(n- 


— -t + ao 

n—rn ' z 


n— m. 


n = m 




«o 


(1- 


m)ait + «2 


a±4> 




,n = 


a\x + ao 


9{t) 




— 2a\(j)/m 


4 

m = 3 


,n = 


a^x 2 + a\x + «o 


g(t) 




-|(2a 2 x + ai)4> 



To apply the classical method to equation Hl.l|) . we consider the one-parameter Lie group of 
infinitesimal transformations in (a;, t, <p) given by 

x* = x + e£(x, t, 4>) + 0(e 2 ), 

t* = t + er(x,t,(j)) + 0(e 2 ), (2.1) 
4>* = 4> + eii(x,t,(l)) + 0(e 2 ), 

where s is the group parameter. Then one requires that this transformation leaves invariant the 
set 

S A = {(f>(x,t) : A = 0} 

of solutions of 1)1. 1J) . This yields an over determined, linear system of equations for the infinite- 
simals £(x,t,cft), T(x,t,4>) and r](x,t,cp). The associated Lie algebra is realised by vector fields 
of the form 

d d d 

v = £0, t,<f>)-Q^ + t(x, t, 0)^ + r](x, t,4>)-Q^- (2.2) 

Having determined the infinitesimals, the symmetry variables are found by solving the charac- 
teristic equation 

dx dt d(f) 



£(x,t,</>) T{x,t,4>) r){x,t,4>y 

which is equivalent to solving the invariant surface condition 

tp = £(x, t, 4>)<t> x + t(x, t, 4>)(pt - r]{x, t, 4>) = 0. 

The set Sa is invariant under the transformation 1)2.1)) provided that pr( 3 )v(A)| a=o = where 
p r (3) v j s ^he third prolongation of the vector field (|2.2|) . which is given explicitly in terms 
of £, r and 7] (cf. |25|). This procedure yields a system of 28 determining equations, a linear 
homogeneous system of equations, for the infinitesimals £(x,t,cj>), r(x,t,(p) and rj(x,t,4>) as 
given in Tabled these were calculated using the MACSYMA package symmgrp .max 



3 Painleve ODE test on travelling waves 

In this section, we use the travelling wave ansatz 4>(x,t) = if)(z), where z = x — ct to reduce 
equation 1)1. 1|) to the ordinary differential equation 

cV>V" + c(n - 3m)^V" - cm(n -m- l)(ip') 3 + n?/> m+ V - cV> m_n+ V = 0, (3.1) 
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Table 2. Results of the Painleve test for Re (p) < 0. 



m, n 



P 



a 



m > 1, m + n > 2 
m > 1, n > 1 



m > 1, n = 1, c / 1 



i 



1— m 
2 

1— m 
-2 



m— n 
1 



m > n, n < 1, m + 2 
m > 1, n = c = 1 
m > 2, n = c = 1 



2-m 
2 

1— m 



arbitrary 



2c 



m— 1 



l/(m-l) 



2(2-m-n) 
(m— n) 2 

arbitrary 
arbitrary 

2c 

(m-l)(l-c) 



l/(m— n) 



l/(m-l) 



1 m— 1+n 
J ' m-1 



2n 
m— 1 



2(2-n-m) 



m—n ' m—n 



'■> m-1 

1 —2771 



m-2 



m— 1 



where (•)' denotes differentiation with respect to z. We then apply the Painleve ODE test as 
described in [2], but in order to do this, we must rewrite 1)3. 1J) in the form 

with the function F analytic in z and rational in its other arguments. If necessary, a transfor- 
mation must be made before applying the test. 

3.1 Integer cases 

We will first start by restricting our investigation to m and n being integers so that the criteria 
on F are satisfied automatically. We now follow [2] and substitute the ansatz 



ip ~ a(z — zo) p , a^O, 

into ()3.1|) . with zq being an arbitrary constant. This gives 

ca 3 p\p(l - m) - l][p(l -m + n)-2](z- z ) 3p " 3 

+ npa m+2 {z - z ) p(m+2) " 1 - cpa m - n+3 (z - ab )P("»-»+3)-i = 0. 



(3.2) 



(3.3) 



It is initially assumed that this is valid in the neighbourhood of a movable singularity and so 
Re (p) < 0. We examine ()3.3)) for all possible dominant balances and for each case, we substitute 



%j)~a{z- z Q ) p + (3(z - Zq) 



P+r 



(3.4) 



into the simplified equation obtained by retaining only the dominant terms. Keeping just the 
linear terms in f3 leads to a cubic equation for r and the solutions of this are known as the 
resonances of the system. If p is not an integer and if the ansatz ()3.2)1 is asymptotic as z — > zq, 
then the leading order behaviour is that of a branch point. Similarly, if the resonances are 
not integers, this also indicates the presence of a branch point. If either of these situations 
occur, then the equation is not of Painleve type and the case can be rejected. The details of 
this procedure are given in Tabled listing p, a and the resonances (apart from r = — 1, which 
always occurs), along with the conditions on m and n for the balance to be dominant. 

We now re-examine ()3.3)l but considering Re (p) > so that ip itself is finite at z = zq 
but one of its derivatives may have a singularity. The results for the dominant balances and 
corresponding values of p and a are presented in Table 01 

Considering the results in Tables [21 and |31 we find that there are ten possibilities for integer 
solutions m and n which give rise to dominant balances with only integer values of p and these 
are listed in Tabic 21 
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Table 3. Results of the Painleve test for Re (p) > 0. 



m, n 



a 



m < 1, m + n<2 

1 + n — m > 0, n > 

m < 1, < 1, ?t, ^ 

777 < 77 , 77 > 1, 777 + 77 ^ 2 

777 < 0, 77 = 

777<1, 77 = C=1 

777 <2, 77 = C=1 

777 < 1, 77 = 1, C 7^ 1 



1 



1— m 
2 



1+n— m 
2 

I— m 
2 

n—m 

-2 
m 

1 

1— m 
2 

2- m 
2 

1—m 



arbitrary 
arbitrary 

-2c 
1— m 



l/(m-l) 



2(2-m-n) 
(m—n) 2 

2(2-m)l 1 / m 
5 — 

m 

arbitrary 
arbitrary 

2e 

(l-m)(l-c) 



1 / (m—n) 



ri 1—m—n 
U ' 1 



0, 



771 

2(m+n-l) 
1+n— m 

2, ^ 

' 1— m 
-2 2(n+m-2) 



n—m' n—m 

2_ 2(2-m) 
m ' m 



o. - 
o, 



-2m 



m-2 



l/(m-l) 



J > l-m 



The resonances in these cases are all integers and they correspond to the points in the solution 
expansion at which arbitrary constants can be introduced. The resonance r = — 1 arises from 
the arbitrariness of zq, whilst if a is undetermined, then r = is another resonance. For each 
pair (p, a), we carry out the solution expansion on the full equation in integral powers of z — zq, 
from (z — zq) p up as far as (z — zo) p+Tmax , where r max is the largest positive resonance. If 
an inconsistency arises, this suggests that logarithmic terms are required in the expansion and 
hence the ordinary differential equation fails the Painleve test. A repeated resonance indicates 
a logarithmic branch point with arbitrary coefficient and this also fails the test 0. Negative 
resonances are not well-understood and if they occur for p < 0, then formally the ordinary 
differential equation fails the test. In a number of the cases, r = —2 occurs in conjunction with 
p = 2, but this then corresponds to a Taylor series beginning with a constant term i.e. p = 0, 
which was not considered and for which there are no inconsistencies in the expansion. Thus, 
these particular negative resonances are spurious and can be ignored. 

There are four integer cases which pass the test, namely m = 0, n = —1; m = 0, n = —2; 
777 = 7i = c = 1 and 777 = 2, 77 = c= 1. However, the last two only pass for one specific value of 
the wave speed c and therefore will not be completely integrable. 



3.2 Rational cases 

When 777 and 77 are not integers, we must first transform the third order differential equation 

so that it is in the correct form for applying the Painleve ODE test. Using the transformation 

i/j = w l , (3.5) 

where I is an appropriately chosen positive integer, we can obtain a suitable ordinary differential 
equation for w, where w'" is a rational function of w, w' and w" . We note that if ift is given 
by (|3.4|) . then the above transformation gives the equivalent expansion for w as 

w^ ai (z-z Q ) pi +p 1 (z-z ) Pl+r \ 

where 



Pi = p/U r i = r - 



(3.6) 
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Table 4. 


Intej 


;er solutions of the Painleve test. 


m n 


V 


a 


r 


D —9 


i 

X 




1, 4 




1 


arbitrary 


0, 3 




2 


-l/2c 


-2, 4 


-1 


-2 


6 


2, 6 




1 


arbitrary 


0, 2 




2 


— 1 /Or 


z, z 


1 


1 


arhitrarv 


0, 




2 


(c-l)/2c, c± 1 


-2, -2 


1 1 
1 — 1 


— i 


_1 L 


1 9 


1 


-2 


2 


2, 2 


1 1 


2 


arbitrary 


0, 2 


1 2 


1 


arbitrary 


0, 2 




2 


-1/2 


-2, 2 


2 1 


-2 


2c/(l -c), c / 1 


-2, 2 




-1 


arbitrary 


0, 2 


2 2 


_2 


2c 


-2, 4 




-1 


arbitrary 


0, 3 




2 


arbitrary 


0, 6 


2 3 


-2 


2c 


-2, 6 




-1 


arbitrary 


0, 4 




1 


arbitrary 


0, 4 



Thus the resonances are left unaltered by this procedure, but the power p is changed by a factor 
of /. We can now use the dominant balances and resonances found for the integer values in 
Tables |U and [3 bearing in mind ()3.6|) . 

We also note that there are no balances in the list which are appropriate for m + n = 2, with 
m > |. This suggests that a simple power is not the correct choice of ansatz and that we need 
to incorporate logarithmic terms at leading order. In fact, we have 

^ ~ a(z - z ) 1/(1 - m) [ln(z - ^)] 1/(2(1_m)) , 

with 

a 2(m-l) = 1 



2(1 - m) 

Thus this case fails the Painleve test due to the logarithms. 

It is now helpful to draw an mn-plot to find which regions of the plane will fail the Painleve 
test. For negative values of p, failure occurs where we have negative resonances and these occur 
when 

m > 1, n > 1, 

n < 1, m + n > 2, 

m > 2, n = 1. 
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Figure 1. The region of the mn-plane which needs to be considered for rational values of m and n. 

For positive p, we can eliminate the cases where p + r < 0, namely 

n > 0, m + n < 0, 

m < n, n > 1, m + n < 1, 

m < 0, n = 0, 

m < —1, n = 1. 

We note that if p/l is to be an integer, then p must be an integer and we can use this to rule 
out further regions in the mn-plane. Thus we have failure when —1 < p < 0, which leads to 

m — n > 2, n<l, m + n^2. 

Similarly for positive values of p, we have failure if < p < 1, and this gives us the two extra 
conditions 

m < 0, m + n < 2, 

m — n < — 1, n>0. 

The remaining region is shown in Fig.Q an d includes the boundaries, the shaded interior and 
the line segment n = 1, 1 < m < 2. This can then be examined in more detail to find the 
rational points which pass the Painleve test. Those which were not examined in the previous 
subsection are found to be 



m 



m 



m 



m 



m 



l 

2 " 
1 

2' 

2 

3' 

3 

2> 

i 

3' 



n = 
n = 0, 
n = 0, 
n = 1 = c 
n = 1 = c 



Again, the last two can be discounted as potential candidates for being completely integrablc 
because they only hold for one value of the wave speed c. 
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3.3 Solitary wave solutions 



In looking for solitary wave solutions of the magma equations, it is helpful to integrate equa- 
tion (|3,1|) twice to yield 



2ip 



l+m 



,1+m— n 



1 



2 ,1+ra-n 2 Alb V 

+ ~ — + -r- + 2Bi/} 2m , 



m 



n 



1 



m 



(3.7) 



for m/l,m+)i/2 and m + n ^ 1, where A and .B are arbitrary constants. Because the scaling 
on equation gave a background level of ip = 1, we restrict our interest to this situation. 

Nakayama and Mason examined the existence of rarefactive solitary waves for which 
ij) > 1 and proved that such waves exist for n > 1 but do not exist if < n < 1. Their proof 
can easily be extended to show that there are no such rarefactive solitary wave solutions for 
n < 0. Nakayama and Mason [22] also looked at the existence of compressive solitary waves 
when < ip < 1 and showed that these were only possible when m = n > 1. However, 
their boundary conditions were unnecessarily strict, permitting algebraic waves but excluding 
exponential waves. 

Here, we give some exact solutions which are possible when the right hand side of ()3.7I) is 
a quartic in ijj. 



(i) m = 0, n = —2. 

?/> = ! + 



3\/2(c + 2) 



^27 - (c + 5) 2 cosh J^(z-zo) -V2(l + 2c) 



This solution describes new exponential compressive waves which exist for —8 < c < —2. 
(ii) m = 0, n = — 1. 

3(c + l) 



V> = 1 



ccosh 



+ c 



This is also a new compressive solitary wave solution when — 3 < c < — 1 and this is 
illustrated in Fig. |2 with c = —2. 

(iii) m = 2, n = 2. 



i> = i 



tp = i 



3(c-2) 



±(1 — c) cosh 
12 

12 + (z-2 ) 2 ' 



2 - z o, 



+ 2c-5 



(3.8a) 



(3.8b) 



Equation ()3.8a|) with the positive sign gives a rarefactive solitary wave which exists for 
2 < c < 4 and agrees with that found by Nakayama and Mason [21 j . However, with the 
minus sign, it represents a new exponential compressive wave solution of amplitude 1 for 
c < and c > 2 and on letting c — ► 2, it gives the algebraic wave (|3,8bj) seen in Nakayama 
and Mason [20j. These waves are shown in Fig. O f° r different values of the wave speed c. 

We now look at other rational values of m and n. By using the substitution (|3.5|) to transform 
equation (|3.7|) into an equation for w of the form 



(wf = f(w), 
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Figure 2. A compressive solitary wave for m = 0, n = — 1 when c = — 2. 



and looking for values of m and n which give f(w) as a quartic polynomial, we can find further 
solutions of this form. It turns out that only I = 2 is possible for real solutions, and the equation 
for w is 



2m 



+ 



W 



2+2m-2n 



+ 



Aw 



2m— 2n 



B 



4m-2 



2c(l — m) 2(2 — m — n) 2c(l — m — n) 

The values of m for which both the first and last terms on the right hand side can give a quartic 
in w are m = | and m = | (with m = 1 being excluded as the denominator of the first term 
would vanish). Choosing appropriate values for A and B so that there is a double root of the 
right hand side at w = 1, we can find solutions for it; and hence i/j. 



(iv) m 



2' 



y, with c < 



2' 



1 + 



1 + 2c 



V 



-^c cosh 



l±2c( 

2c \ Z 



This is a compressive wave which is oscillatory for c < —2 and represents a new solution. 
(See Fig. E|) 



(v) m 



2 1 



n 



0. 



IP 



6 



cosh(z — zq) + 1 

This new solution is oscillatory in nature, has both compressive and rarefactive regions 
and the shape is independent of velocity. (See Fig. [SJ) 



(vi) m 



2- 



n 



1. 







cosh 



6 

-{z - Zq) 



+ 5 
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Figure 3. Solitary waves for m = 2, n = 2, which are exponential for c = 
and algebraic when c = 2 ( ). 



-2 ( ) and c = 3 ( ) 



This is a compressive wave with amplitude 1 for c < or c > 1. 

2 ' " 2 • 

\ 2 



(vii) m = |, n - 



1 + 



2c -3 



2 

12 



cosh 



2c-3 
2c 



- zo) 



+ 2-c 



(z - z ) 2 + 



(3.9a) 
(3.9b) 



Equation (|3.9a() represents rarefactive waves with the positive root if 3 > c > | (found by 
Nakayama and Mason whereas it gives new compressive oscillatory waves for c > | for 
the negative root. As c — ► | for the compressive wave, we find the algebraic solution ()3.9b|) 
which is given in Nakayama and Mason |2(Jj . 

4 Other scaling reductions 

For the five cases which pass the Painleve ODE test on the travelling wave reduction, we now 
look at the test on the other symmetry reductions found in Section [21 

(i) m = 0, n = —2. Using the reduction, 

z = xt 1/2 , cp(x,t) =t 1/2 w(z) 
leads to the equation 

zww'" - 2zw'w" - zw 3 w' + 3ww" - 4(w') 2 + Aw' - w 4 = 0. (4.1) 
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Figure 4. An oscillatory solitary wave for m = — A, n 



i when c 



(ii) m = 0, n = —1. Here the appropriate reduction is 

z = xt 1/3 , 4>(x,t) = t 2/3 w(z), 
leading to the equation 

zww'" - zw'w" - zw 2 w' + 4ww" - 3{w') 2 + 3w' - 2w 3 = 0. 

(iii) m = i,n = — i. In this case, we use 

z = xt 1/2 , <f>(x,t) =tw 2 (z) 
to obtain 

zww'" — zw'w" — zw 3 w' + 3ww" — 2{w') 2 + w' — w 4 = 0. 

(iv) m = 2, n = 0. The correct reduction for this case is 

z = xh(t), 4>(x,t) = x~ 4 w 2 (z) 
where h is an arbitrary function and this leads to the equation 
z 2 w"' - 2zw" - ww' + 2w' = 0. 

(v) m = |, n = 0. Lastly, we have the reduction 

z = xh(t), (f)(x, t) = x~ 3 w 3 (z) 
yielding 



(4.2) 



(4.3) 



(4.4) 



z 2 w"' - w 2 w' = 0. 



(4.5) 
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Figure 5. A solitary wave with compressive and rarefactive regions for m = |, n = 0. 
Table 5. Results of the Painleve test for Re (p) < 0. 



m n p a r 






-2 


-1 


±V2 


1, 


4 





-1 


-2 


6 


2, 


6 


1 

2 


1 

2 


-1 


±2 


2, 


4 


1 

2 





-2 


12 z 2 


4. 


6 


2 
3 





-1 


±V6z 


3, 


4 



As in the travelling wave reductions, we again use the ansatz 

ip ~ a(z — zo) p , a^O, 

to find the dominant balances for Re (p) < and the resulting resonances, with r = — 1 always 
occurring due to the arbitrariness of zq. The results are shown in Table [S] and give integer values 
for p and r. 

However, for cases (|4.1j) . (|4.4j) and (|4.5|) . the solution expansion on the full equation as far 
as powers of (z — zo) p+rmax gives inconsistencies, indicating that logarithmic terms are required 
in the expansion. Therefore, these cases fail the Painleve test. This leads to only cases (14.21) 
and 1)4 .3|l remaining, which both satisfy the condition m = n + 1 found by Harris . 

5 Painleve PDE test 

In this section, we consider the two remaining cases and relate them to known completely 
integrable equations. 
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(i) m = 0, n 



1. A slightly more general equation is considered, namely 



(5.1) 



which reduces to the magma equation for m = 0, n = — 1 upon setting fj, = 1. 

It is straightforward to apply the Painleve PDE test, as described by Weiss, Tabor and 
Carnevale jHUj, to Q5.1[l and it is found to pass the test for any values of \i. We note that (|5.1|) 
can be written in the form 



Then making the transformation <j> = v x and integrating once gives 
Vxxt ~ v x v t - /X = 0, 

where we have set the function of integration to zero, without loss of generality. A simple 
transformation 

v(x, t) = —u(x, t) — x — t 

then leads to the shallow water wave (SWWI) equation 

u xxxt 4" U x U x t + UtU xx — u x t — u xx = 0, (5-2) 

which was considered by Clarkson and Mansfield This is also equivalent to the Hirota- 
Satsuma equation 



where d~ 1 (f)(x) = f(y) dy, which was discussed by Hirota and Satsuma [T4*| . 

Clarkson and Mansfield jOj have discussed symmetry reductions of (|5.2jl obtainable using 
the classical Lie method and the nonclassical method of [I]. They gave a catalogue of classical 
and nonclassical symmetry reductions and exact solutions of (|5.2|) . Of particular interest are a 
plethora of solutions of (|5.2|) possessing a rich variety of qualitative behaviours. These arise as 
nonclassical symmetry reductions and all of which look like a two-soliton solution as t — > — oo, 
yet are radically different as t — > oo. These results have important implications with regard 
to numerical analysis and suggest that solving (|5.2j) numerically could pose some fundamental 
difficulties. An exponentially small change in the initial data yields a fundamentally different 
solution as t — > oo. How can any numerical scheme in current use cope with such behaviour? 
(See jH] for further details.) 

The Lax pair for the Hirota-Satsuma equation 1)5. 3 J) is the third scattering order problem [S] 



[{</> xt - //) /4>] x -<k = 0- 



Vxxt + VVt - Vxd x 



(Vt) ~Vx-Vt = 0, 



(5.3) 



Ipxxx - (l + \v) tp x = Xtp, 

3\^ t = [l-^imMxx- 



r\ti> x - 



(5.4a) 
(5.4b) 



We remark that (|5.4a|) is similar to the scattering problem 



ipxxx + 3(1 + 6u)tp x + I u x - iV3d x 1 



(ut) ^ = Xtp, 



which is the scattering problem for the Boussinesq equation 



tlxxxx ~\~ 3 (it )xx U X x — l^tti 

and has been comprehensively studied by Deift, Tomei and Trubowitz [5]. 
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(ii) m = 21 n = ~\- Similarly, we consider a slightly more general equation 

uu xxt - u x u xt - u 3 u t + \iu x = 0, (5.5) 

where setting 4>(x, t) = u 2 and \x = \ yields the magma equation for m = ^, n = — ^. 

Again, the Painleve PDE test due to Weiss, Tabor and Carnevale |3flj can be applied and ([5.5)1 
is found to pass for any /x. Clarkson, Mansfield and Milne [7] showed that (|5.5|) arises as 
a symmetry reduction of the (2 + l)-dimensional Sine-Gordon system 

( &xt \ ( ®yt \ ®yRxt - @xRyt _ q 



ysinGy^ \sin0yy sin 2 9 

(Rxt\ (Ryt\ @y@xt - 9 
\sin@ ) x \sin9/ y sin 2 9 



derived by Konopelchenko and Rogers |16U17j . The special case of (|5.5(l with fj, = is equivalent 
to the sine-Gordon equation 

u xt = sinn, (5.6) 

which is one of the fundamental soliton equations solvable by inverse scattering using the AKNS 
method pQ. 

Clarkson, Mansfield and Milne show that a Lax pair associated with 1)5.5)1 is given by 
i\) x = {— iA<T3 — ^iu<T2} "0, (5.7a) 
= < JTT! ""3 + tt^i " ^72^2 \i> = 0. (5.7b) 



4Xu 4A 8A 2 

where <ti, <T2 and cr^ are the Pauli spin matrices given by 

01\ / -i \ /10 



4X1 = 1 i o > CT2= l i o h CT3 = V -1 



Equations 1)5.7)1 are compatible, i.e. tp xt = tp tx provided that u(x,t) satisfies 1)5.5)1 . We remark 
that the Lax pair ()5.7)l reduces to that for the sine-Gordon equation 1)5 .6)1 if /i = p. Further 
we note that the spectral problem ()5.7a)l is the standard AKNS spectral problem whilst if \i ^ 0, 
then l)5.7b)l involves powers of both A -1 and A~ 2 . 



6 Discussion 

We have performed a comprehensive Painleve analysis of the generalized magma equation 

fa + (<p n + mr-^Mx - 4> n ~ m 4>t x ) x = o, 

and found that there are only two pairs of values of the indices m and n for which the equation is 
completely integrable. The case m = 0, n = — 1 is related to the Hirota-Satsuma equation and 
was studied by Clarkson and Mansfield [Bj. The Hirota bilinear form, multi-soliton solutions 
and the Lax pair are known for this equation. The situation when m=^,n = — ^ is less clear 
cut. It is related by a simple transformation to a special case of the real, generalized, pumped 
Maxwell-Bloch system examined by Clarkson, Mansfield and Milne [7|. It has an non-isospectral 
Lax Pair and a single soliton solution, but the Hirota bilinear form and multi-soliton solutions 
are as yet unknown. 
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